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FLIGHT-DETERMINED AERODYNAMIC STABILITY AND CONTROL DERIVATIVES 


OF THE M2-F2 LIFTING BODY VEHICLE AT SUBSONIC SPEEDS 

Robert W. Kempel and Ronald C. Thompson 
Flight Research Center 


SUMMARY 


The longitudinal and lateral-directional stability and control derivatives determined 
from flight tests of the M2-F2 lifting body research vehicle are presented for a range 
of Mach numbers from 0. 41 to 0. 64 and for angles of attack from -2. 3° to IS. 8°. The 
derivatives were determined by using a high-speed repetitive analog-matching technique 
and approximate equations to analyze recorded time histories of vehicle response to 
random control inputs and pulses. The flight-determined derivatives were compared 
with wind-tunnel predictions for representative Mach number and angle-of-attack con- 
ditions. Using these derivatives, the vehicle modal-response characteristics were 
calculated. 

Comparisons of flight and wind-tunnel results indicated that the flight-determined 
longitudinal aerodynamic static stability, lower-flap effectiveness, and pitch damp- 
ing were generally higher than predicted. The normal-force-curve slope as well as 
the lateral-directional aerodynamic stability and control effectiveness generally agreed 
with wind-tunnel predictions. In general, the lateral-directional-damping derivatives 
were in good agreement with predictions. 


INTRODUCTION 


The first manned flight tests of lifting body reentry vehicle configurations were 
made at the NASA Flight Research Center, Edwards, Calif. , under the joint sponsor- 
ship of NASA and the U. S. Air Force. The test program was designed to investigate 
performance, stability and control, and flying-quality requirements of this new class 
of vehicle throughout the subsonic, transonic, and low supersonic flight regions. Test- 
ing was started with the lightweight M2-F1 vehicle; the results of these tests are re- 
ported in reference 1. Subsequently, to expand the limited performance envelope of 
the basic M2-F1, a heavier and more sophisticated vehicle, the M2-F2, was designed 
and built. The M2-F2 was designed for launch from a B-52 carrier aircraft and was 
capable of either a glide flight to landing or an accelerated flight to higher altitudes and 
speeds utilizing a modified LR-11 rocket engine, followed by a glide flight to landing. 
The first glide flight of the M2-F2 was made on July 12, 1966; in the following 8 months, 
16 glide flights were completed. On the last flight, before the powered phase of the 
program, the vehicle was damaged extensively in a gear-up landing. 


Even though the powered phase of the program was not conducted, the initial glide 



tests served many important functions. Prior to the flight program, predictions of the 
flight characteristics of lifting bodies were, by necessity, based almost exclusively on 
wind-tunnel data. The flight program did, therefore, make it possible to compare the 
extensive small-scale and limited full-scale wind-tunnel results with the actual flight 
data in order to establish some measure of prediction accuracy. In addition, a hybrid 
computer program was mechanized for the determination of flight derivatives based on 
a curve-fitting technique which used the data obtained from the flight program. The 
flight results also made possible an assessment of the sensitivity of handling qualities 
to typical discrepancies between predicted and actual flight derivative values. 

In this report the more significant longitudinal and lateral-directional aerodynamic 
stability and control derivatives obtained from flight are compared with small- and 
full-scale wind-tunnel results. The modal-response characteristics calculated from 
the flight-determined derivatives are included. 

General stability and control flight results for the M2-F2 vehicle are reported in 
reference 2, and performance, lift, and drag results are reported in reference 3. 

ta> 


SYMBOLS 


Data are presented as standard NASA coefficients of forces and moments which are 
referred to the vehicle body axes fixed with respect to the vehicle center of gravity. 
Positive directions are: X, forward; Y, to the right; and Z, down. Positive direc- 
tions of the forces, moments, and angular displacements and velocities are in accord 
with the right-hand rule. 


Physical quantities are given in the International System of Units (SI) and paren- 
thetically in U. S. Customary Units. The measurements were taken in U. S. Customary 
Units. Factors relating the two systems are presented in reference 4. 
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VEHICLE DESCRIPTION 


General 

The M2-F2 was a single-place, wingless, research vehicle (figs. 1 and 2). The 
design was basically a blunt, 0. 227 radian (13°) half cone with a flat upper surface and a 
boattailed afterbody. Two vertical fins were mounted high on the afterbody and equipped 
with trailing-edge rudders. The upper trailing edge of the afterbody was provided with 
two flaps, and the lower afterbody trailing edge was provided with a single-section 
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Figure L M2-F2. 





Figure 2 . Three-view drawing of the M2-F2 vehicle. Dimensions in meters (feet) unless otherwise indicated. 


flap (figs. 2 and 3). The structure incorporated a semimonocoque construction with 



Figure 3. Rear view of the M2-F2. 


two full-depth keels. The physical characteristics of the M2-F2 are listed in tables 
1, 2, and 3. 


Flight Control System 

The flight control system of the M2-F2 was an irreversible, dual, hydraulic system. 
Conventional stick and rudder pedals were provided for roll, pitch, and yaw control, 
and artificial feel was provided about all axes. Cockpit control displacements and 
gradients and corresponding control surface travels are summarized in tables 1 and 2. 

Coarse pitch trim was provided by positioning the upper flap by means of a trim 
wheel in the cockpit. For fine pitch trim, the pilot positioned the lower flap by means 
of a trim switch on the control stick. 

Pitch control was provided through longitudinal displacement of the control stick 
which positioned the lower flap. Roll control was provided through lateral displacement 
of the control stick which differentially positioned the upper flaps. Yaw control was 
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TABLE 1. - PHYSICAL CHARACTERISTICS OF THE M2-F2 VEHICLE 


Body — 

Planform area, meters^ (feet2): 

Actual 

Reference, S 

Longitudinal length, meters (feet): 

Actual 

Reference 

Span, without rudder flare, meters (feet): 

Actual 

Reference, b 


Aspect ratio, 


b£ 

S ’ 


basic vehicle 


14. 9 (160) 

12. 9 (139) 

6.76 (22.2) 

6 . 11 ( 20 . 0 ) 

2. 94 (9. 63) 
2. 91 (9. 54) 

0. 655 


Body-leading-edge sweep, degrees 77 

Lower flap — 

Area, meters^ (feet^) 1.41 (15.23) 

Span, meters (feet) 1. 65 (5. 42) 

Chord, meters (feet) 0. 86 (2. 81) 

Deflection, degrees: 

Pilot's control authority, down 5 to 30 

Pitch stability augmentation system authority ±5 

Upper flaps, two — 

Area, each, meters^ (feet^) 0. 89 (9. 57) 

Span, each, meters (feet) 1.31(4.28) 

Chord, meters (feet) 0.68 (2.23) 

Deflection, each flap, degrees: 

Pitch trim (symmetric travel), up 0 to 35 

Pilot's aileron authority (asymmetric travel) ±5 

Roll stability augmentation system authority 

(asymmetric travel) ±2 l/2 

Vertical stabilizers, two — 

Area, each, meters^ (feet^) 1. 50 (16. 10) 

Height, trailing edge, meters (feet) 1. 16 (3. 79) 

Chord, meters (feet): 

Root 2.24 (7.36) 

Tip 0. 79 (2. 58) 

Leading-edge sweep, degrees 62. 3 

Rudders, two — 

Area, each, meters^ (feet^) 0.49 (5.27) 

Span, each, meters (feet) 1. 28 (4. 20) 

Chord, meters (feet) 0.38 (1.25) 

Deflection, each (outward), degrees: 

Pilot's effective control authority 12 

Yaw stability augmentation system authority 4. 2 
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TABLE 2.- FLIGHT CONTROL SYSTEM OF THE M2-F2 


— 

Surface 
rate 
limits, 
deg/ sec 

25 

30 

22 
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provided through displacement of the rudder pedals. The surfaces operated about a 
5° flared condition as illustrated in figure 2. In response to a pedal input, one surface 
moved out as the other moved in, although the rudder could not move inboard farther 
than the faired or 0° deflection point. 

Because of the high adverse yaw created by aileron deflection, the M2-F2 was 
equipped with an aileron-to-rudder interconnect; that is, a mechanical system which 
deflected the rudder by an amount proportional to the aileron deflection. The desired 
ratio of rudder-to-aileron setting was determined by the pilot and was controlled 
through a wheel in the cockpit. This feature enabled the pilot to make coordinated turns 
by using aileron inputs only. For these tests the interconnect ratio was set at a nominal 
value of approximately -0. 5. 


Stability Augmentation System 

A rate-feedback stability augmentation system (SAS) was incorporated about all axes 
to provide angular rate damping. The SAS mode and gains were controlled by the pilot 
by means of cockpit switches. Roll and yaw channels with high pass filters provided 
the capability of canceling or washing out SAS control surface inputs when the angular 
rates approached a constant. Constant rate rolling and yawing motion was therefore 
not impeded by SAS control surface deflections. 

Weight, Balance, and Mass Distribution 

Precise M2-F2 weight and balance measurements were made whenever vehicle 
weight and mass distribution changed significantly. These measurements then became 
the reference condition used for the digital computer program in the update computation 
of mass distribution characteristics for various loading conditions. This program com- 
puted vehicle weight, three-coordinate center-of-gravity locations, moments of inertia, 
product of inertia, and inclination of the principal axis. Computations were made by 
compiling the weight and location of each item of structure, hardware, instrumentation, 
expendable material, ballast, and miscellaneous equipment in the vehicle with respect 
to the reference axes. It was also necessary to correct the program by introducing a 
bogie weight which brought the computed weight and longitudinal center-of-gravity lo- 
cation into agreement with the actual reference weight and balance measurements. Ta- 
ble 3 presents the computed weight, three -coordinate center-of-gravity positions, and 
moments of inertia of the M2-F2 used in this analysis. 

Before the M2-F2 was completely outfitted for flight, it was decided to determine 
experimentally the yawing moment of inertia, the product of inertia, and the inclination 
of the principal axis at a lightweight condition as a check on the digital computer pro- 
gram. Table 4 compares the experimentally determined weight, inertias, and centers 


TABLE 4. - COMPARISON OF EXPERIMENTAL AND COMPUTED WEIGHT AND MASS 
DISTRIBUTION CHARACTERISTICS OF THE M2-F2 


Quantity 

Experimental 

Calculated 

W, kg (lb) 

Iz, kg- m2 (slug-ft2) 

IXZ> kg- m2 (slug-ft2) 

Inclination of principal axis, deg 

Vertical center of gravity below reference, cm (in. ) 

Horizontal center of gravity, cm (in. ) 

15X9 (4976) 
8080 (5958) 
-638 (-471) 
-5.3 

30. 8 (12. 1) 
341 (134. 1) 

1519 (4976) 
8025 (5925) 
-569 (-420) 
-4.8 

37. 6 (14. 8) 
341 (134. 1) 
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of gravity with computed values. It should be noted that the bogie weight was generally 
approximately 2 percent of the weighed reference value and never exceeded 3 percent. 
On the basis of these results the digital program was used with a high degree of confi- 
dence. 


FLIGHT TESTS 


The M2-F2 was air-launched from a B-52 aircraft at approximately 13, 730 meters 
(45, 000 feet) altitude and a Mach number of 0. 65. All flights were unpowered glide 
flights and were made under strict visual flight rule (VFR) conditions. All landings 
were on Rogers dry lakebed, which is approximately 700 meters (2280 feet) above mean 
sea level. 

All maneuvers were performed at angles of attack between -2. 3 ° and 13. 8 ° and at 
Mach numbers from 0. 41 to 0. 64. They consisted of the conventional pulse-type ma- 
neuvers used for stability and control analysis and large angle-of-attack excursions 
used for performance, lift, and drag analysis. Both lateral-directional and longitudinal 
SAS-off pulse maneuvers were performed for specific identification of aerodynamic 
stability derivatives. SAS-on maneuvers, with controls not held fixed, were also per- 
formed for analysis of aerodynamic stability and control derivatives. For the analysis 
of any single M2-F2 flight response, the time interval was generally from 5 to 8 seconds. 
Short times were used primarily to maintain tolerable limits from^ relative steady-state 
conditions. 


INSTRUMENTATION 


Data were acquired by means of a pulse code modulation (PCM) system. No pro- 
visions were made for recording data onboard the vehicle. All data were sampled 200 
times per second, conditioned, and transmitted to the ground station. On the ground 
the data were recorded on magnetic tape and processed, utilizing appropriate cali- 
brations and corrections. Vehicle translational accelerations and angular rates were 
not corrected for instrument location because this factor was determined to be negligi- 
ble. Airspeed and altitude were measured by using a NASA pitot static tube mounted on 
the nose boom (fig. 1). The static- and total-pressure orifices were, respectively, 

1. 704 meters (5. 59 feet) and 1. 905 meters (6. 25 feet) ahead of the vehicle nose. The. 
angle of attack and angle of sideslip were measured by free-floating vanes on the nose 
boom (fig. 1). The angle-of-attack vane was 1. 20 meters (3. 95 feet) ahead of the nose 
of the vehicle. The angle-of-sideslip vane was approximately the same distance from 
the nose of the vehicle, but oriented 90° to the axis of the angle-of-attack vane. Cor- 
rections for upwash and sidewash were applied to the angle-of-attack and angle-of- 
sideslip data. (See ref. 3. ) 

The basic accuracy of the PCM system was generally within 2 percent of full-scale 
deflection. The accuracies of the various parameters used in this report are presented 
in table 5. 


13 



TABLE 5 .- ACCURACIES OF PERTINENT PARAMETERS 
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METHODS OF ANALYSIS 


An analog-computer-matching analysis procedure was used in this study. This pro- 
cedure is generally applicable to any aerodynamic vehicle, and particularly to motion 
which is separable into longitudinal or lateral-directional components with continuous 
control inputs. In the analog-matching technique, the equations of motion representing 
the mathematical model were programed on an analog computer. The computer was 
"forced" to solve the equations of motion based on the vehicle flight-recorded control 
surface inputs (forcing functions) as a function of time. Each analysis consisted of a 
maneuver at a specific time and time interval for a particular set of flight conditions. 

All flight-recorded variables of interest were stored in a digital computer and displayed 
to a human operator on a cathode ray tube (CRT). The computer was operated at a rate 
faster than real time so that the computer solutions and recorded flight time histories 
could be displayed on the CRT simultaneously as standing wave forms. By manually 
adjusting the aerodynamic stability and control derivatives of the computer-mechanized 
mathematical model, the operator attempted to match the computer-generated response 
with the measured vehicle flight response. This process of curve fitting, or difference 
minimization, allowed an experienced operator to make judgments and tradeoffs on the 
basis of phase and amplitude relationships and to manually adjust the derivatives accord- 
ingly. He could also disregard information which appeared to be incorrect. A match 
was considered to be achieved when a minimum difference was obtained simultaneously 
between all the flight-recorded and computed variables. More detailed information con- 
cerning the process of extracting aerodynamic stability derivatives from flight results 
is included in references 5 to 8. 

The equations of motion used throughout this study were linearized, rigid-body 
equations referenced to vehicle body axes. These equations, which are developed in 
references 9 to 13, are presented in the appendix. The longitudinal calculations were 
generally based on two-degree-of-freedom equations with a provision for linear vari- 
ations of dynamic pressure when conditions dictated. The longitudinal parameters 
matched were q, 0, cu, and a n . Three-degree-of-freedom lateral-directional equa- 
tions were used. Lateral-directional parameters matched were p, <p, 1 3, r, and ay. 

When longitudinal or lateral-directional, pulse-type, SAS-off maneuvers were per- 
formed, simplified approximate analyses were used to corroborate the results of the 
matching process. The simplified equations used in the computation of the approximate 
derivatives are also presented in the appendix. 


RESULTS AND DISCUSSION 


Longitudinal Static Stability and Control 

The operational flight corridor of the M2-F2 is shown in figure 4 in terms of upper 
and lower limits of dynamic pressure. The data points indicate conditions for which 
longitudinal maneuvers were analyzed. 

The least troublesome mode, in terms of derivative extraction, was the longitudinal 
short-period mode. The dampers-off, short-period-pulse, transient response maneu- 
vers matched were determined to be essentially linear, thus validating the use of 
linearized equations in the analysis. A typical pulse-induced, longitudinal, short -period, 
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2394.0N/m z 


launch point 


q =19, 150 N/m z ^X\ 

ILlCl^t x\\ 


Edwards altitude 


Figure 4 L M2-F2 flight corridor , indicating conditions for which longitudinal maneuvers were analyzed. 

dampers-off, flight transient response and the computer-generated, matched time 
history are presented for comparison in figure 5. This is considered to be a relatively 
good match. 







Figure 6 is a semilogarithmic plot of the half-amplitude envelope of the flight- 
measured transient response of figure 5 (At is measured from the time where 6^ re- 
turns to its static position). The exponential decay in the response variable amplitude is 
typical of the response of linear systems. By using the approximate literal expressions 



Figure 6. Amplitude versus time for the example flight-measured transient longitudinal response of figure 5. 
T j j 2 = 1.45 sec; P = 2.73 sec; £u) n - 0.478 rad/sec; |^| =2.26 per sec; co n = 2.31 rad/sec. 


presented in the Longitudinal Equations section of the appendix 
mation from figure 6, the derivatives C m 


ex’ 


'' m 6i 


and 


m. 


in conjunction with infor- 
+ C m ^\ were obtained 


for comparison with the matched derivatives. Table 6 compares the analog-matched 
and approximate derivative solutions for this flight condition. The agreement is 


TABLE 6. - COMPARISON OF AERODYNAMIC DERIVATIVES DETERMINED BY 
ANALOG MATCHING AND BY SIMPLIFIED CALCULATIONS FOR 
ONE LONGITUDINAL PULSE 



Cm a 

(S + Cm o) 

C™ 

Analog matching 
Simplified calculations 

-0. 00169 per deg 
00146 per deg 

-0. 492 per rad 
666 per rad 

-0. 00247 per deg 
00263 per deg 
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considered to be good, thus validating the use of the approximate derivative approach 
when pulse-type response data are available; however, the analog- matched derivatives 
are considered to be more correct because fewer approximations were used in obtaining 
them. 

An example of a typical longitudinal flight maneuver time history with a random- 
appearing input and analog match is presented in figure 7. This is also considered to 
be a relatively good analog match. 




Figure 7. Flight-measured and analog-generated longitudinal random input time history. SAS on; 
V = 156 mjsec (512 ft/sec); q = 3480 N/m 2 (73 Ib/ft 2 ); M = 0.53. 


19 



Table 7 summarizes the longitudinal flight data analyzed and the results obtained 


TABLE 7.- FLIGHT-MEASURED LONGITUDINAL AERODYNAMIC DERIVATIVES AND CALCULATED RESPONSE CHARACTERISTICS 


q , 

or, 

deg 

V, 

V 

M 

« u . 

deg 

* 1 . 

deg 

c m a . 
per deg 

^ tn q + c m^. 
per rad 

Cm v 

per deg 

c N , 

INqj 

per deg 

^sp> 

rad/sec 

rad/sec 

N/m 2 

lb / ft 2 

m/sec 

ft/sec 

m 

ft 

5837 

121 9 

5.8 

ESI 

441 

6,308 

20,697 

0. 43 

-11.3 

19.4 

-0.00121 

-0. 547 

-0. 00243 

0. 0232 

2. 089 

0. 518 

5841 

122. 0 

5.3 

BS8S 

438 

6,004 

19,697 

.42 

-11.4 

19.5 

-. 00146 

-. 657 

-. 00265 

.0249 

2.228 

. 606 

7129 

148 9 

2.9 

H89 

566 

9,584 

31,443 

. 60 

-11.7 

21.2 

-.00182 

-.475 


.0267 

2.773 

.464 

7824 

163.4 

1.2 

156.0 

512 

6,246 

20,492 

.49 

-11.4 

21.2 

-. 00129 

-.479 

IT 1 9 

.0293 

2.740 

. 613 

8714 

182.0 

2.2 

170. 1 

558 

6,846 

22,459 

.54 

-11 3 

22.0 

-. 00145 

-. 507 


.0277 

2.442 

. 591 

5937 

124. 0 

4.7 

182.3 

598 

11,219 

36,809 

.62 

-11 4 

20.0 

-.00169 

-. 492 

-. 00247 

.0294 

2.436 

389 

4209 

87.9 

10 7 

163. 7 

537 

12,009 

39,400 

.56 

-11.3 

16.4 

-.00151 

-.446 

-. 00237 

.0283 

1. 925 

.281 

8090 

169. 0 

2.0 

170.4 

559 

7,507 

24,630 

55 

-8. 5 

16. 2 

-. 00174 

- 400 

-. 00257 

.0316 

2. 882 

507 

8044 

168 0 

2.0 

157.9 

518 

6,244 

20,486 

.50 

-11.4 

21.5 

-.00164 

-. 654 

-. 00238 

.0312 

2. 777 

. 747 

9030 

188. 6 

2. 1 

164.3 

539 

5,857 

19,215 

.52 

-14 2 

25.0 

00162 

-. 579 

-. 00219 

0304 

2. 816 

484 

8556 

178.7 

1.2 

174. 6 

573 

7,634 

25,045 

.56 

-11 4 

20.7 

-.00157 

- 415 

lEBliIlWiM 

.0249 

2. 932 

.734 

3754 

78.4 

13. 0 

162. 5 

533 

12,654 

41,517 

55 

-12 4 

11.5 

-. 00098 

-.400 

-. 00183 

.0236 

1.384 

.204 

7800 

162.9 

1.4 

150.0 

492 

5,586 

18,327 

.47 

-14 2 

21.5 

-. 00129 

-.560 

-. 00179 

. 0234 

2.284 

. 568 

8259 

172.5 

4.2 

160 3 

526 

6,149 

20, 174 

51 

-14.3 

19. 0 

-.00132 

-.400 

-. 00173 

.0258 

2. 389 

.465 

3495 

73. 0 

13 4 

163. 7 

537 

13,137 

43, 101 

.56 

-14 8 

15 0 

- 00125 

-. 475 

-. 00188 


1.510 

.216 

3481 

72.7 

13. 5 

156. 1 

512 

12,471 

40,916 

.53 

-14.8 

14. 9 

-. 00108 

- 400 

-. 00198 


1. 399 

. 199 


Figure 8 is a comparison of M2-F2 wind-tunnel data and flight-determined static longi- 
tudinal derivatives. The wind-tunnel derivatives presented in this report were obtained 
from reference 14 and from unpublished data. Ail wind-tunnel data were adjusted to 
correspond to the flight vehicle center of gravity, which was approximately 54 percent 
of the reference vehicle length of 6. 1 meters (20 feet). Full- and small-scale wind- 
tunnel data were measured about a moment reference center located at 55 percent of the 
reference length. 




Figure 8. Comparison of M2-F2 wind-tunnel and flight-determined static longitudinal derivatives. 


20 




































Values of the flight-determined derivative at the low angles of attack were 

generally higher than predicted; at medium and higher angles of attack, the agreement 
was good. The flight-determined lower-flap effectiveness, C m ^ , was generally higher 

than comparable wind-tunnel data except at high angles of attack where agreement was 
good. The flight-determined normal-force-curve slope, , was also generally in 

good agreement with tunnel predictions. The scatter band in the flight data of figure 8, 
although large in some of the data, would have been reduced if more data had been avail- 
able and if better-conditioned flight maneuvers could have been performed. However, 
the results of this study are believed to be representative of a flight-test program of 
this type. 

The longitudinal modal-response characteristics were calculated by using equations 
similar to the longitudinal equations in the appendix. The equations were solved by 
using a digital computer program in which the matched derivatives were the inputs. 
These data are presented in table 7. 

Lateral-Directional Static Stability and Control 

The operational flight corridor and related conditions for which lateral-directional 
flight maneuvers were analyzed are presented in figure 9. A typical example of a 


2(M03 



Figure 9. M2-F2 flight corridor indicating conditions for which lateral-directional maneuvers were analyzed . 
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flight-measured, lateral-directional maneuver with a continuous random input and its 
computer-matched time history is presented in figure 10. Figure 11 is a time history 
of a dampers-on, aileron-doublet pulse showing the flight response and computer match. 


10 

6a, deg o 


-10 

40 

p, deg/sec 0 






Flight data 

Analog match 


Figure 10. Flight-measured and analog-generated lateral-directional random input time history with varying aileron- 

to-rudder interconnect ratio. SAS on; a = 3.2 °; V = 196 m/sec (642 ft/sec); q = 6910 N/m^ ( 144.4 lb/ft^ ); 

M = 0.63. 
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Flight data 







Figure 1L Flight-measured and analog-generated lateral-directional aileron-doublet input time history . SAS on; 
a = 3.1°; V = 164 m/sec (537 ft/sec); q = 7280 N/m 2 (152 lb/ ft 2 ); M = 0.53. 


Figure 12 presents a typical Dutch roll mode, dampers-off, transient response flight 
time history and the resulting computer match. Each of the computer-matched time 
histories of figures 10, 11, and 12 shows acceptable agreement with its respective 
flight response. 
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. 4 

6 r , deg 

2 

0 12 3 4 5 

t, sec 

Figure 12. Flight-measured and analog-generated lateral-directional Dutch roll mode time history of a pulse maneuver . 
SAS off; a = 5.8°; V= 179 mjsec (587 ft /sec); q = 6040 N/nt 2 (126 lb, /ft 2 ); M = 0.61. 



Figure 13 is a semilogarithmic plot of the half-amplitude envelope of the flight- 
measured transient response of figure 12 (At is measured from the time when 6 r re- 


turns to its static position). This plot was used in the simplified analysis and shows 
good linearity. The equations for the simplified approximations are contained in the 
Lateral-Directional Equations section of the appendix. The derivatives determined 
from the simplified approximations are compared in table 8 with the results from analog 
matching for the time history of figure 12. For this example, the directional-stability 
derivative, C n ^, and the rudder-effectiveness derivative, C n g , determined from 

approximate equations show good agreement with the analog-matched values; however, 
the effective-dihedral derivative, , is approximately 12 percent lower than the 

analog-matched value. It is interesting to note that the approximation for the effective- 
dihedral derivative is highly dependent on the measurement of the ratio of roll rate to 


sideslip angle. For the M2-F2 this ratio was generally between 3 and 15. The peak-to- 
peak amplitude of the sideslip angle measured in flight seldom exceeded 2°, and the effect 
of sidewash on the sideslip vane was not well defined; therefore, sideslip angle might be 


suspected as a source of possible error which could affect the 



value determined. 


A distinct advantage of the computer-matching technique is suggested in the fore- 
going discussion. If a particular parameter is needed but not required to accomplish a 
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Figure 13. Amplitude versus time for the example flight-measured lateral-directional Dutch roll mode response of 

[pi | r I 

figure 12. Tjj 2 = 2.04 sec; = 0.34 rad /sec; u> n ^ = 4.80 rad/sec; P = 1.31 sec; |-| =14.2; |^| =3.3. 


TABLE 8. - COMPARISON OF AERODYNAMIC DERIVATIVES 
DETERMINED BY ANALOG MATCHING AND BY 
SIMPLIFIED CALCULATIONS FOR A 
LATERAL-DIRECTIONAL PULSE 



C^, per deg 

Cry per deg 

c n 6 > Per de S 
°r 

Analog matching 
Simplified calculation 

-0. 00796 
00700 

0. 00608 
. 00628 

-0. 00224 
00205 
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computer match, and that parameter is suspected of being in error or is possibly not 
available, the operator would adjust the credibility of that parameter accordingly in the 
matching process. That is, the operator would weight each parameter of the time 
history on the basis of the credibility of each and would make proper tradeoffs between 
the parameters in the time history which were known to be correct and those that were 
in error. The best compromise would be one that would result in the best set of deri- 
vatives for any particular case. It should be noted, however, that a best match does 
not necessarily represent the best set of derivatives. 


Table 9 summarizes the lateral-directional responses for which analog matches 
were attempted and the results obtained. Figures 14(a) and 14(b) are plots of lateral- 
directional wind-tunnel and flight-determined static-stability and control derivatives. 
The wind-tunnel derivatives were obtained from reference 14 and unpublished data. 


The flight values of the dihedral derivative, 


C 1 ^, the aileron-effectiveness derivative, 


Ci . , and the rolling-moment-due-to-rudder-de flection derivative, C , are in good 
°a 

agreement with wind-tunnel results (fig. 14(a)). Flight values of the directional-stability 
derivative, C n ^, the rudder-effectiveness derivative, C 


n 6. 


and the yawing- moment- 


due- to-aileron-de flection derivative, 



are also in good agreement with wind- 


tunnel predictions (fig. 14(b)). It was not generally possible to detect differences in the 
aerodynamic derivatives over the Mach number range (0. 42 to 0. 64) of these tests, 
although the possibility that Mach effects were present should not be precluded. 


The lateral-directional modal-response characteristics were calculated by using 
equations similar to those in the appendix. The equations were solved by using a 
digital computer program in which the matched derivatives were used as inputs. These 
data are presented in table 9. One interesting point indicated by this study was the 
relatively high ratio of bank-angle-to-sideslip-angle amplitude for the Dutch roll mode 
at low angles of attack due to the high value of the effective-dihedral derivative. Also 
of interest is the predicted existence of a coupled roll-spiral mode. Actual verification 
of this mode was difficult because of its strong dependency on the roll-damping deriva- 
tive, C£p, and the cross-derivative, Cnp. It should be noted that the response 

characteristics presented in table 9 are for the unaugmented (no dampers) vehicle and 
that the addition of rate damping would usually break this mode into the more common 
roll and spiral modes. 


Longitudinal and Lateral-Directional Damping and Cross -Derivatives 

Generally, the most difficult derivatives to predict are the cross-derivatives and 
the damping derivatives. Predicted values for the M2-F2 were obtained from a com- 
bination of theoretical estimates and unpublished wind-tunnel data. These predictions 
are believed to represent the best average values over the Mach and angle-of-attack 
range of interest. Also, because the analog-matching process is less sensitive to 
variations of the cross- and damping derivatives than to variations of any other deriva- 
tives, the damping derivatives are more difficult to obtain from flight responses. For 
example, in the analog-matching process it was possible to vary the predicted damping 
derivatives from 25 to 50 percent without seriously affecting a match. The cross- 
derivatives could have been changed from 50 to 200 percent without affecting a match. 
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(a) Rolling-moment derivatives . 

02 1— i 



O Flight data 

Wind tunnel 

Full scale (ref 14), M = 0.2 5 

Small scale, M =0.6 



(b) Yawing-moment derivatives. 

Figure 14. Comparison of M2-F2 flight and wind-tunnel lateral-directional rolling- and yawing-moment derivatives. 
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Values of the pitch-damping derivative, + C m ^, are shown in figure 15. The 

average flight-determined pitch-damping derivative was generally from 25 to 50 percent 
higher than predicted. This fact was corroborated by pilot comments concerning the 
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Figure 15. Comparison of M2-F2 flight-measured and predicted damping and cross-derivatives. 

short-period-mode damping from flight pulses as compared with a flight simulation of 
the M2-F2 in which the predicted values were used. 

The lateral-directional-damping derivatives and cross-derivatives are also shown 
in figure 15. All data used in the matching process are included, although in some 
instances the data were not changed from the predicted value because, as previously 
stated, moderate variations from these values had little effect on the individual match. 
In general, the roll-damping derivative, C; , and the yaw -damping derivative, C n „, 

P r 

which were determined from flight data are in good agreement with predictions. 


Side- Force Characteristics 

The relative insensitivity of the transverse accelerometer used in the flight program 
precluded the determination of any reliable value of the side-force derivative, Cy^. 

As for the cross- and damping derivatives, the side-force derivative was difficult to 
obtain from the analog-matching process. As previously mentioned, the peak-to-peak 
amplitude of the sideslip angle seldom exceeded 2°. At a dynamic pressure of 
14,400 N/m^ (300 lb/ft*), 1° of sideslip would produce approximately 0. lg transverse 
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acceleration. Figure 16 compares the wind-tunnel and flight-determined side-force 
derivatives. Only the data points which were considered valid are presented. 



Figure 16. Comparison of side-force derivatives obtained from flight and wind-tunnel results. 


CONCLUSIONS 


Longitudinal and lateral-directional aerodynamic stability and control derivatives 
for the M2-F2 vehicle were extracted by using an analog computer program which 
matched the flight responses of the vehicle with corresponding computer responses. 
Analysis of the data from the first 16 glide flights revealed the following: 


1. At low angles of attack the longitudinal static-stability derivative, C mQ; , was 


generally higher than predicted by wind-tunnel results. The lower-flap-effectiveness 


derivative , C 


m 6 


l 


was also higher than predicted except at high angles of attack 


where agreement was good. The normal-force-curve slope, was in relatively 

good agreement with the predicted values. 


2. The lateral-directional static-stability and control-effectiveness derivatives 
were in good agreement with the wind-tunnel-predicted values. 


3. In general, the lateral-directional cross- and damping derivatives were in good 
agreement with the predicted values. However, the longitudinal pitch-damping deriva- 
tive was 25 to 50 percent higher than predicted. 

Flight Research Center, 

National Aeronautics and Space Administration, 

Edwards, Calif., July 21, 1971. 
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APPENDIX 


EQUATIONS USED FOR DERIVATIVE DETERMINATION 


Longitudinal Equations 

Equations used in the mechanization of the analog-computer-matching process were 
as follows: 

Normal-force equation 

V(a - q) = Z w w + Zgj 6 1 - (g sin 0 O ) 9 
Pitching-moment equation 

q = M w w + Mqq + M^w + 

The increment of normal acceleration as measured by a body-axes -oriented acceler- 
ometer is expressed as 

Aa n = | [<* - q + g (sin 0 O ) oj 

Simplified approximations used in determining longitudinal static derivatives and 
damping derivatives are 

(■% E p) 2 “ZwM,-VMw 


% sp“”sp”" (Zw+ M 0 + VM i> 

„ W ^ a n 
Im ^ — 

o; qS A a 


C 



I X /Aq \ 

qS cW/ 


Lateral-Directional Equations 

Equations used in the mechanization of the analog-computer-matching process were 
as follows: 

Rolling-moment equation 

ixz 

P = j-“ r + Lp/3 + L p p + L r r + L 6 __6 a + L 6 r 5 r 
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APPENDIX 


Yawing-moment equation 



+ N^/J + NpP + N r r + N 6a 6 a 4- N 6r 6 r 


Side-force equation 


f - P“o - r + (|) V + V s + y 6 a «a + *8 r «r 


The increment of transverse acceleration as measured by a body-axes -oriented 
accelerometer is expressed as 




Simplified approximations used in determining lateral-directional effective dihedral 
and directional-stability derivatives and the rudder-effectiveness derivative are 
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